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Abstract 

Tlie determination of the Fermi coupling constant, Gf, is examined in the light of re- 
cently calculated 2-loop QED corrections and planned experiments to measure the muon 
lifetime to a level below 1 ppm. The methods used in the calculation of the QED cor- 
rections are described in detail. Sources of the dominant theoretical and experimental 
uncertainties are identified. Finally the incorporation of into analyses using the full 
electroweak Standard Model is discussed. 
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1 Introduction 



Any renormalizable theory only becomes predictive once it has been supphed with a sufficient 
number of experimental inputs so as to fix the free parameters that appear in its Lagrangian. 
In order for the subsequent theoretical predictions to be as accurate as possible, these inputs 
are chosen from the experimentally best-measured quantities available. In the case of the 
Standard Model of electroweak interactions these are the electromagnetic coupling constant, 
a, the Fermi coupling constant, Gp, and the mass of the boson, Mz- Their recognized 
best values, along with their absolute and relative errors that they represent are [0,|2| 

a = 1/(137.0359895 ± 0.0000061) (0.045 ppm) 

Gf = (1.16639 ± 0.00002) x 10"^ GeV"^ (17ppm) 
Mz = 91.1867 ± 0.0021 GeV (23 ppm) 

The first of these, a, is measured at momentum q = and when used in the analysis high 
energy data is afflicted with a 'hadronic uncertainty' entering already at the 1-loop level 
and arising because it must be 'run up' from low energy crossing, on the way, the hadronic 
resonance region. 

Most processes, such as e~^e~ annihilation, that allow Mz to be accurately determined, 
are also mediated by the photon and thus contain a pair of characteristic energy scales. 
Experiments performed at low-energy therefore do not provide useful insight into physics at 
the high scale, Mz- 

The Fermi coupling constant is obtained from the muon lifetime, r^, via a calculation in 
the Fermi model, in which the weak interactions are represented by a contact interaction. It 
suffers from neither of the disadvantages described above. There can be no hadronic effects 
at 1-loop in this model and even in the full Standard Model 1-loop hadronic effects in muon 
decay are suppressed by a factor m-j/M^ where is a light fermion mass. Charged current 
processes are only mediated by the W boson whose effects are unobscured no matter at 
what energy experiments are performed. The measurement of the muon lifetime can thus be 
considered a high-energy experiment that is actually performed at a low scale This point 
is highlighted by the fact that the next generation of muon lifetime measurements will be 
sensitive to W propagator effects. 

In the mid-80's, just before the turn on of LEP, a CERN report concluded that the error 
on Mz would be ±50 MeV or 550 ppm and that "A factor of 2-3 improvement can be reached 
with a determined effort" Q. It was thus generally thought that Mz would be the input 
parameter limiting the accuracy with which theoretical predictions could be made. As is 
often case, due to a variety of unforeseen effects, the uncertainty turned out to have been 
significantly overestimated and it now approaches that Gp. The lesson that should be 
learned from this is that it is extremely difflcult to predict, even in the relatively short term, 
the accuracy to which fundamental parameters will be determined and it is important that 
these be extracted to the limits that the current theoretical and experimental technology 
allows. A great effort and expense was brought to bear in order to reduce the uncertainty on 
Mz to the level given in the table above and whereas it is unlikely to diminish significantly 
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in the near future a muon collider offers the prospect of a reduction by a factor of 10 

The Fermi coupling constant is closely related to the p-parameter that was introduced by 
Ross and Veltman as a way of interrogating the mass generation sector of the theory 

-l + 5p (1.1) 



" M| cos2 Ow 

The detailed understanding of the nature of the Higgs mechanism will come about through 
the scrutiny of systems and processes that are most sensitive to the Yukawa sector of the 
theory. The Fermi coupling constant obtained from the muon lifetime provides such a probe 
but stands alone for the tremendous accuracy with which in can be measured. For this reason 
a considerable amount of theoretical effort, reviewed in section |8|, has been devoted to the 
study of higher order weak corrections to the muon lifetime. It can be shown that the 2-loop 
QED corrections to the muon lifetime affect the predicted value of the Higgs mass at the 
percent level which is precisely the level needed to probe the detailed structure of the Higgs 
system by its radiative corrections. 

In this paper we look closely at the theoretical and experimental considerations necessary 
for the extraction of Gp from to an accuracy of 1 ppm or better as is anticipated for new 
experiments planned at the Brookhaven National Laboratory, the Paul Scherrer Institute 
and the Rutherford-Appleton Laboratory. What is known about the relationship between 
the muon lifetime and the Fermi coupling constant is reviewed in section ^ In section ^ 
inadequacies that appear at higher orders with the usual definition oi Gp are pointed out 
and the relative merits of the one adopted here are discussed. In section Q it is shown how 
to perform the renormalization of the electromagnetic coupling constant in such a way as to 
incorporate the large logarithms, ln(mg/m^). Section H describes the methods used for the 
calculation of the 2-loop quantum electro dynamic (QED) corrections to the muon lifetime 
in the Fermi model. The sources of the dominant theoretical and experimental uncertainties 
are examined in sections |6| and |^ respectively. The incorporation of the value oi Gp obtained 
here into analyses using the full electroweak Standard Model is discussed in section ^. 

Appendix ^ contains expressions for certain loop integrals for which exact analytic forms 
are known in dimensional regularization. Appendix ^ gives the results for the individual 
Feynman diagrams that occur in the calculation of the 2-loop QED corrections to muon 
lifetime and Appendix]^ gives an expression for the electron spectrum to 0{a) keeping the 
full electron mass dependence. Finally Appendix ^ examines the branching ratio for the 
process /i~ — > e~t'^^'ee^e~. 

2 The Muon Lifetime and The Fermi Couphng Constant 

In the Minkowskian metric in which time-like momenta squared are positive, the Lagrangian, 
relevant for the calculation of the muon lifetime in the Fermi theory is 

= Cqeb + ^QCB + (2.1) 
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Here jO-q^j^ is the usual bare Lagrangian of Quantum Electrodynamics (QED), 

^°QED = E - - 7 M - - E Qfi^hpi^'K- (2-2) 



The sum is over all fermion species, /, with wavefunction, -0^-, mass, mj, and electric charge, 
Qf. Ap \s the photon field. The superscript " indicates bare, as opposed to renormalized, 
quantities. ^Cq^j-, is the bare Quantum Chromodynamic (QCD) Lagrangian responsible for 
strong interactions. The Fermi contact interaction that mediates muon decay is 



-w 



-2V2Gf [^°,7a7lV'°] . [4°7A7L<] (2.3) 



in which ipe, tpUfj, and ip^^^ are the wavefunctions for the muon, the electron and their asso- 
ciated neutrinos respectively, and 'Jl denotes the usual Dirac left-hand projection operator. 
For the present purposes the Fermi coupling constant, Gp, goes unrenormalized. 

To leading order in Gi;' and all orders in a the formula obtained for the muon lifetime, 
T^, by means of the Cp takes the general form 



r^ = ro(i + Ag). (2.4) 



where 



To = — 2.5 

" 1927r^ ^ ^ 

and Aq encapsulates the higher order QED corrections and can be expressed as power series 
expansion in the renormalized electromagnetic coupling constant ar = e,^/(47r). 



Aq = J2^1^'^ (2-6) 



i=0 



in which the index i gives the power of that appears in Ag^*). 

Although Cp is not renormalizable, the Ag^*-* can be shown to be finite 0. This fol- 
lows from the fact that when Fierz rearrangement is used to rewrite C\y in so-called charge 
retention order. 



Cw ^ -2V2Gf [V^°7a7l<] • 7a7lV': 



the currents remain purely V— A in form. By contrast, in the case of neutron decay, the anal- 
ogous transformation generates scalar and pseudo-scalar terms and the following arguments 
break down. The radiative corrections in that case are not finite. Considering the vector 
part ipej^iipfj. of this effective fi-e current, one sees that after fermion mass renormalization 
is performed the remaining divergences are independent of the masses and thus cancel, as 
for the case of pure QED. The QED corrections to the axial vector part may be shown to 
also be finite by noting that the transformations V'e 75V'e and mg —rUe leave >Cqed and 



4 



£qcd invariant but exchange 'tpelx^P^i ^ V'eTATsV'/i- Thus the radiative corrections to the 
axial-vector contribution to the muon decay matrix element can be obtained from those of 
the vector contribution by setting nie —me. 

This argument was used by Roos and Sirlin P to show that terms odd in nie would cancel 
between vector and axial-vector contributions in the expression for the differential decay rate. 
They then went on to show, by direct examination of known analytic form of the differen- 
tial decay rate |^], that the phase-space integration could not generate terms linear in me 
and thus that the leading electron mass corrections at 1-loop are O (a(mg/m^) ln(me/m^)) , 
O (a(mg/m^)) and higher but they did not exclude the possibility of O [a{m^/m^)) terms 
which, in fact, do occur. 

It is possible to show that Aq cannot contain any terms that are odd in me at any order in 
a. The total decay rate may be calculated directly as the imaginary part of muon self-energy 



diagrams and was done in Ref. [10|. In charge retention order terms in the numerators that 
are odd in the electron mass lead to a helicity flip along the internal electron line that causes 
the purely left-handed V— A vertices to annihilate. However the integrals themselves generate 
non-analytic terms such as ml ^Jva^ = |mep which may appear to be odd in me if the absolute 
value is dropped. 

The above considerations are true in any regularization scheme and have the importance 
consequence that in the limit m-g — > only the radiative corrections to the vector pieces in 
Hw/ need to be calculated as they are equal to those to the axial vector part. This avoids 
entirely the problems associated with 75 when dimensional regularization is used. 

Dropping the electron neutrino mass and keeping only the leading term in that of the 
muon neutrino gives 

Ag(°) = -8x-12x2lnx + 8x^-x^-8y + C'(xy), a: = — 7, V = (2.7) 

that comes from phase space integrations. Notice that the coefficient of the leading terms in 
the electron and muon neutrino masses are identical which again follows from the fact that 
their wavefunctions may be interchanged by Fierz rearrangement of £w 
The first order corrections to Ag are 

Ag(i) = (^^) - 3C(2) - (34 + 12 In x)x + 96C(2) xi + 0(x In^ x) j (2.8) 

where C is the Riemann zeta function and (^(2) = vr^/G. The leading electron mass-independent 
term was calculated Kinoshita and Sirlin |11] and it was the observation that the result is 
well-behaved in the limit mg — > that ultimately lead to the discovery of the Kinoshita-Lee- 
Nauenberg (KLN) theorem ||l^. An exact expression for the full electron mass dependence 
in Ag(^) has been given by Nir |13|. Starting from the expressions for the 1-loop QED 



corrections to the electron spectrum given by Behrends et al. ||9[ and, after taking into account 
the correction given in Appendix C of Ref. |jll[, we obtain complete agreement with Eq.(12) 



of Ref. |13|. As a spinoff we have considerably simplified the expression for the electron 



spectrum of Ref. [Q. The result appears in Appendix |y. An analogous expression has been 
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given by Czarnecki et al. |14| for the QCD corrections to 6 ^ crvr- We have also checked 
the electron mass-dependent terms in Eq. (|2.8| ) by performing a large momentum expansion 
of the imaginary part of propagator-type diagrams. 

The second order corrections to /S.q have been presented recently |1^, 15] for mg = 0. The 
result, ignoring the effects of tau loops, is 

V TT / V 5184 27 ^ ' 36 ^ ^ 8 ^ ^ 

+ 53C(2) In 2 -(0.042 ±0.002)) (2.9) 



where ^(3) = 1.20206... and C(4) = 7r^/90. The numerical constant is the hadronic contri- 
bution with a conservative estimate of its error. The contribution from tau loops has been 



shown to be very small [15| as anticipated from the decoupling theorem 

(2) 



TT 



Hl = -[ — ] 0.00058 (2.10) 



and will be neglected unless otherwise stated. 

The methods used to obtain the result ( |2.9| ) will be described in some detail in section |^. 
A key feature is that the calculation was performed by means of cutting relations ||l^ applied 
to the 25 1-particle irreducible (IPI) diagrams that appear in Fig.|^. In these diagrams the 
thick line represents a muon with the external legs being on-shell. The thin lines represent 
the electron, electron neutrino or muon neutrino all of which are taken to be massless. The 
wiggly line represents the photon. Cuts passing only through massless internal lines give 
non-zero contribution and all others vanish. In particular, since the external, leg is on-shell, 
any cut through a muon line is identically zero. The imaginary part of the diagrams of Figjl] 
then generates precisely those combinations of the products of amplitudes that appear in the 
calculation of the 2-loop QED corrections to muon decay. Moreover, when the calculation 
is performed in this way, the intricate cancellation of infrared (IR) divergences, that would 
occur between these products of amplitudes, coming from distinct cuts applied to a given 
diagram, is largely taken care of automatically. 

The result for /S.q^'^^ in Eq. (2.9) is composed of several independent pieces. The part 



coming from the purely photonic diagrams of Fig.|l| that contain no closed charged-fermion 
loop and related 1-particle reducible (IPR) external leg corrections is 

. /ar\2 /11047 1030 , , 223 , , 67,, , , .\ 

— 3.55877. (2.12) 



vr 



This was calculated in a general gauge and was found to be both gauge invariant and 
ultraviolet (UV) finite. 

Feynman diagrams containing an electron loop or e'^e" pair in the final state are obtained 
from cuts to diagrams D6-D8 of Figjl|. Addition of the related IPR external leg corrections 
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Figure 1: 4-loop diagrams whose cuts give contributions to the muon decay rate. The 
diagrams are grouped according to the main integration topologies A-G (see Figs. The 
results for the individual diagrams are given in Appendix ^ 
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and diagrams containing coupling constant counterterm insertions gives 



Or 



vr 



2 



— 3.22034 (2.14) 



The value obtained in Eq.( 2.14 ) is consistent with a numerical study presented in Ref. [0] in 
the context of semi-leptonic decays of heavy quarks, Q — > XgePg- 

The contribution from diagrams containing muon loops, diagrams C3-C5 of Fig.|l] can 
be shown to be identical in the on-shell and MS renormalization schemes for a 't Hooft 
mass /i = (see section^). These plus related IPR external leg corrections and diagrams 
containing coupling constant counterterm insertions and give 

—) 0.0364333. (2.16) 

vr / 

This was calculated both by means of dispersion relations in the on-shell scheme [|15| and 
by the perturbative methods described in section |5| in the MS renormalization scheme. Full 
agreement was found. 



The hadronic contribution to (|2.9| ) is given in Ref. |15|. Notice that, because it was 
calculated using dispersion relations, it naturally involves a subtraction of the photon vacuum 
polarization at q = which corresponds to the on-shell scheme being adopted for electric 



charge renormalization. The same is true for the stated result for tau loops of Eq. ( 2.10| ). On 



the other hand the photonic and electronic contributions, ( |2.11 ) and ( 2.13| ), were calculated 



perturbatively assuming, where necessary, the MS renormalization scheme. It will be shown 
in section ^ that the two sets of results can be combined in a consistent manner by setting 
renormalized electromagnetic coupling constant, Ur = ae(?7^^)=l/135.90. In doing so all 
corrections up to ©(a^), O (a^ ln(mg/m^)) and O (a* ln*^"'^(mg/m^)) for all i > 2 are taken 
into account. 

Fermion mass renormalization is always performed in the on-shell renormalization scheme 
in which the renormalized mass of a stable fermion is set equal to its pole mass. 

Aq^"^^ of Eq. ( |2.9D , when combined with the tiny contribution of tau loops gives Eq. ( 2.10| ), 

Ag(2) = To (^^^^^ ' (6.700 ± 0.002). (2.17) 



When this is used in Eq. ( |2.4D along with the current best value for the muon lifetime, 
= (2.19703 it 0.00004) ^s M, it gives a new value for the Fermi coupling constant 



Gf = (1.16637 ± 0.00001) x lO^'^ GeV"^ (9 ppm) 

This represents a reduction in the overall error on of about a factor of 2 and a downward 
shift in the central value of twice the experimental uncertainty. The error is now entirely 
experimental. 
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3 The Definition of the Fermi CoupUng Constant 



At high precision it is essential to have a clear and unambiguous definition for the Fermi 
coupling constant, Gp- For the present purposes Gp will taken to be defined by Eq.(E^) 



.-1 



1927r3 



1 + Ag) 



(3.1) 



with Ag being calculated using the Fermi Theory Lagrangian Cp oi Eq.(| 



Gp to arbitrary accuracy in terms of the physically observable quantities, r^, a, 



This expresses 
and 

me- In this way the infrared structure, which involves both real bremsstrahlung and virtual 
corrections, is entirely relegated to a single constant Ag and does not have to be dealt with 
when Gf is eventually used in the calculation of other physical constants such as the W 
boson mass, M^. Thus the, presumably well- understood, contribution of QED to the muon 
lifetime is eliminated from Gp which is left with an enriched contribution from weak physics. 
Gp IS sometimes defined via the formula 




1 + 



2tt 



25 
T 



vr 



(3.2) 



where 



and 



F[x) 



l-8x 



12x^ lnx + 8x^ -x^ 



a{m^) ^ = a ^ 



In 



m. 



+ 



67r 



136. 



In as far as this can be derived from Cp \\, provides an adequate definition but it contains a 
number of features that become out of place at higher orders. 

The function coming from phase space integration, does not factorize in this way at 
higher orders. 

The factor (l + (3m^)/(5M^)) is the effect of the W boson propagator and is not gen- 
erated by the Cp. It is more naturally included along with the weak corrections in Ar as 
is described in section |8|. Its presence is an historical artifact of an attempt to reconcile the 
observed muon and neutron beta decay rates with universality of weak interactions before 
the advent of the Cabibbo angle |7,|l8|. In that scenario a light W boson with a mass slightly 



heavier than the kaon was needed so as to forbid the unobserved decay mode, K 



The expression for a{m^) used in the above definition contains a term l/(67r) that comes 
from W boson loops in the photon self-energy. As these do not come from Cp they should be 
omitted so as not to risk double counting when is used as input in electroweak calculations. 

It is sometimes suggested that Eq.( |3.2D should be viewed as an exact definition oi Gp. 
This has significant drawbacks and hard to justify now that the 2- loop QED corrections are 
available [10, lH|. It is already the case that these 2-loop QED corrections are larger than 
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the current experimental error and are roughly an order of magnitude greater than those 
anticipated in the next generation of experiments (see section 0). If Eq. (3^) is taken to be 



exact then these 2-loop corrections will ultimately have to be incorporated in formulas that 
relate Gp to other electroweak observables such as the mass of the W boson, Mw. This is 
extremely inconvenient as one would prefer to deal with quantities from which QED has been 
eliminated as far as possible. 

4 Renormalization of the Electromagnetic Coupling Constant 

In this section it will be shown how to set up the MS renormalization scheme so as to 
obtain the leading logarithmic corrections to the muon lifetime of O (a^ ln(m^/mg)) and of 
O (a* ln*~^(m^/mg)) for all i > 0. It is also shown how to incorporate hadronic contributions 
|jl5[, calculated in the on-shell renormalization scheme, in a consistent manner. 

In any self-consistent calculation using the Lagrangian (^]^), the quantity a that appears 
in the Ag^*) is the renormalized electromagnetic coupling constant, to be denoted a,., in 
whatever renormalization scheme has been chosen. It is generally convenient to adopt the on- 
shell renormalization scheme for the fermion masses which will be assumed throughout unless 
otherwise stated. No such restriction is placed on the coupling constant renormalization. 

For problems containing widely disparate scales, renormalization of the coupling constant 
in the MS scheme is to be preferred as it automatically absorbs the dominant logarithmic cor- 
rections into ar at the outset and avoids the need for resummation of large logarithms coming 
from higher-order contributions as is required when the on-shell renormalization scheme is 
adopted. 

In QED the numerical value of is by obtained solving the equation 

a = ^ . (4.1) 

i-n!;^(o) 



where a on the right-hand side of Eq.(4.1) is the experimentally-measured quantity, a = 
1/137.0359895(61) |^]. Il'^^{0) is the photon vacuum polarization function which itself may 
be written as expansion in a^, 

oo 

n;,(o) = J^n;«(o) (4.2) 

1=1 

and each term receives contributions from all fermion species. The hat, ^ , indicates that 
counterterm contributions in the chosen renormalization scheme have been included. In the 
on-shell renormalization scheme the counterterms are adjusted so that 11^^(0) = and it 
follows from Eq.(^) that the renormalized coupling constant in this scheme satisfies, = 
Q!OS = ct. In the MS renormalization scheme the counterterms are chosen to contain only 
divergent pieces plus certain uninteresting constants. In particular 1-loop counterterms are 
just proportional to A = 1/e — 7^; -|- ln47r -|- 0(e) where D = 4 — 2e is the dimensionality of 
spacetime and = 0.57721566... is Euler's constant. In this scheme an appropriate choice 
for the 't Hooft mass is /i = and we write = a{m^). 
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Defining U'-iji} (0) = aJ.P'^*^ and solving Eq. (|4l|) for a,, yields 

ar = a- a^P^^^ - a^P^^^ + O(a^) (4.3) 

The contribution of leptons to ll''lj{0) can be exactly calculated in perturbation theory 
using dimensional regularization []l9|p0| . The 1-loop contribution is 

fiS'(O) = -J,4,)-fr(2-f)X:Q?(^)*" + 2^|l (4,4) 

where Se^^^ is the leptonic contribution to the 1-loop electromagnetic charge counterterm. 
At this order the MS and on-shell renormalization schemes differ only in the the finite parts 
of their coupling constant counterterms, and ScqI respectively. Notice that the con- 

tribution from muon loops to (^e^ and Scq^ is identical for jj. = and thus their overall 
contribution to n^y^,^'*(0) vanishes in both schemes. 



The 2-loop correction in a general renormalization scheme can be obtained from Ref. |21| 
taking into account the difference in measure used to define dimensional regularization 

^ ^ '-";-^-_",;--> (4.).-^r (3 - f ) r (. - f ) E (f " ' 

where 6m^^^ is the 1-loop lepton mass counterterm and Sef^^ is the leptonic part of the 2- 
loop electromagnetic charge counterterm. Adopting the on-shell renormalization scheme for 
fermion masses gives 



and substituting this in Eq.(|4.6|) yields 



"^V(O) = 4^ D{D-3){D-5) ^^"^ 



X 



a, 

4tt' 



2 / ^2 ^K\ 



^ Qf (a - In ^ + ^) + 2^ + 0{D - 4) (4.9) 
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whose leading logarithms are in agreement with the well-known result of Jost and Luttinger 
|22|. The leading logarithms of C(a^) were obtained by Rosner |23|. 



The contributions to muon lifetime from hadron, muon and tau loops have been calculated 



using dispersion relations [15|. This method naturally generates a subtraction at momentum 
q = which corresponds to on-shell renormalization of the coupling constant. The calculation 



of Ref. 1 15 1 was performed by effectively replacing the photon propagators 



+ ie 



9fio 



+ ie 



+ ie 



(4.10) 



The on-shell results of Ref. |15| are converted to MS by adding a correction obtained by 
replacing the photon propagator by 



^ -|- it 



-|- ie 



)- f<5eW 
e V MS 



6e 



(1) 
OS 



q'^ + ie 



(4.11) 



where — (5eQg is a finite constant that includes the effects of loops of all fermion species 



,(1) 



except the electron. As usual e is a positive infinitesimal. When ( 4.11 ) is applied to the 
calculation of the muon lifetime the terms on the right hand side proportional to g^^i/ cancel 
between diagrams hence its effect amounts to an overall multiplicative factor for the photon 
propagator and leads to a correction to r„ of 



AT 



OS-^MS 



6e 



(1) 

OS 



This can be accounted for by a redefinition of the renormalized coupling constant, of 
Eq.( |4.3| ), which has the effect of eliminating the contributions of all fermion loops except 
those of the electron. It follows that the results of Ref. ||l^ can be used directly in the 
calculation of the muon lifetime provided that the MS value a(/x) is replaced by 



a(/i) — > Oei^i) = a + 



1 a \^ I5a-^ 
Sit Att'^ J ml IGvr^ 



(4.12) 



The KLN theorem [12| requires that the only electron mass singularities contributing the 
muon lifetime are those that appear in (|4T^ ) which arise as a consequence of a's having been 
defined at the exceptional momentum, q^ = 0. Moreover, careful examination of Eq.( [4.1[ ) 
reveals that writing 



a 



+ 



a 1^ '^^J. 



a 

4^ 



2 

In-^ 

mi 



(4.13) 



correctly resums logarithms of the form a" ln"'~^(m^/mg) for all n > and incorporates those 
of 0(0^ ln(m^/mg)) without generating spurious higher-order logarithms as there would be 
if the second term on the right hand side of Eq. ( [4.13|) where moved into the denominator of 
the first. The constant term of 0{a^) has been dropped since terms of this order have not 
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been accounted for in the present calculation. Terms that are not logarithmically enhanced 
are correct up to 0(0^). 

For the extraction Gp in the Fermi theory the appropriate choice is /i = which gives 

ae(m^) = 1/135.90 = 0.0073582 

however when Gp is used in the analysis of electroweak data obtained at the weak scale 
/U = Mz is generally a more convenient choice. 

The expression for ae{m^) could also be obtained be adopting a hybrid renormalization 
scheme in which electron loops were renormalized in MS and all other fermions loops in the 
on-shell scheme. This is allowable in QED since the individual fermions are not connected by 
any gauge symmetry and provides a simple, if unconventional, way of deriving Eq.( 4.13| ). 



5 The Calculation of the 2-loop QED Contributions 

In this section the evaluation of Feynman diagrams that appear in the calculation of Aq^'^^ 
Eqs. (^.6|), (|2.9| ) is discussed. For the calculation of the 25 muon decay diagrams of Fig.| it 
is necessary to deal with only 7 basic 4-loop topologies labeled A-G (see Figs.|3|-|5|). These 
7 basic topologies are assumed to have scalar propagators raised to arbitrary (positive and 
negative) integer powers. 

The integrals are evaluated on the mass shell = m'i using dimensional regularization 



1 19,^ for both the UV and IR divergences. Within dimensional regularization a few integrals 
are known in a simple closed form for arbitrary powers of the propagators as a ratio of Euler 
r functions. These are listed in Appendix |A[ 

The imaginary parts must be obtained for the 4-loop topologies. Those 4-loop integrals 
that, on general grounds, have no imaginary part will be discarded along the way. This 
includes 4-loop vacuum bubble integrals and 4-loop integrals with a through-going on-shell 
line, as these have no cuts that can generate an imaginary part. 

To reduce the burden of analytic integration as much as possible we use the well-known 
method of integration-by-parts in dimensional regularization to express all integrals of a 
given topology, having arbitrary integer powers of the propagators, in terms of a small set 
of primitive integrals. Using integration-by-parts has the advantage of requiring the explicit 
calculation of only a fixed number of primitive integrals instead of hundreds of scalar integrals 
that generally appear in a multiloop calculation. On-shell integration-by-parts relations were 
used at the two- loop level for the first time in Ref . (2^ (but see also ) and at the three-loop 
level in Ref. (2|. 

A given set of primitive integrals is not unique as one can eliminate one set in favor of 
another set of integrals, as long as they are independent with respect to the integration-by- 
parts identities. Therefore only the total number of primitive integrals is a fixed quantity for 
a given topology, and what is used as a primitive integral is a matter of convenience. 

Integration-by-parts relations follow from calculating the derivative in the identity 

d^P^ r/(p,^i,-- - ,ln)] = (5.1) 
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t 0,0,0 





tm,0,0 



' m,m,0 



Figure 2: Examples of triangle diagrams with different mass configurations for which the 
triangle relation Eq. (|5.3|) holds. The lines are numbered according to Eq. (|5.2|). 



within dimensional regularization. Here p is any momentum that is integrated over, k £ 
{p, h, - ■ ■ , In} and / is a scalar function that may depend in addition to p on a set of momenta 
k. 

An important example of an integration-by-parts identity is the triangle relation for one- 
loop 3-point integrals [p3] . Let us define for masses a, b, c 



7'a,fe,c(ai, "2, as, "4, 05, ae) = y d^P ^ 



1 



(p2)"i [{p - kY - a2]«2 [{p - ly - 62]«3 



(5.2) 



Here and below +ie in the propagators is not written explicitly to keep expressions compact 
(e.g. is understood to be l/(p^ + ie)). By taking the derivative in the integral 



d^p 



d 



one obtains the triangle relation 



{D — 2qi — 02 — 03) — 022^(1" 



a33+(l- -4-) 



T'a,6,c(ai,a2,a3,a4,a5,a6) = 

(5.3) 



where the notation of Ref. |2J] has been adopted. Thus the operator 2"*" raises the power 02 
in T by one, 4~ lowers 04 in T by one, etc.. By repeated use of this relation on integrals T 
with positive integer powers , ae one is able to lower at least one of the powers ai, 04 

or 05 to zero. Instead of applying the triangle relations recursively, the result can be obtained 
directly and expressed as a 3-fold nested sum [^8|. Note that the configuration of masses in 
Eq. (^]^) is not completely general, and if all propagators carry different masses additional 
terms appear in the triangle relation, leading to a less useful identity. However, when a 
triangle with the mass configuration of T appears as a subgraph of a multiloop integral, one 
can use the triangle relation to simplify the multiloop integrals considerably. Some subgraphs 
relevant to the present calculation are given in Fig]2| 
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Figure 3: Integration topologies A^^^ and A. The lines are numbered according to Eg. ( [5. 5]) 
and Eq. (5.7) respectively 



To express multiloop integrals of a certain topology in terms of primitive integrals, one 
generally needs to consider all integration-by-parts identities obtained from Eq. ( |5.1| ) and 
solve them. To solve these relations there are two possible ways to proceed^. Integration-by- 
parts relations may be explicitly applied to a large set integrals with specific integer powers 
{• • • , —2, — 1, 0, 1, 2, • • • } and solved the resulting set of linear equations (see, for example, Ref. 
p7|). Alternatively the relations may be considered as symbolic equations with unspecified 
powers of the propagators and combined into a set of recurrence relations that, after repeated 



application, express integrals in terms of primitives |24,25|. 

Both approaches can have advantages in particular circumstances. The second approach 
is more commonly used as one may encounter integrals with many different powers in appli- 
cations that involve expansions of diagrams. For the present work the recursive approach was 



used and a recurrence scheme was implemented in a FORM [3C] program. As some steps in 
the recurrence relations involve hundreds of terms we limit ourselves to the main conclusions 
for each of the basic 4-loop topologies and give a list of the primitive integrals were chosen. 
The calculation of these primitive integrals is discussed in the next section. 

For the 3-loop topology A^^^ of Fig.^, it is possible to lower the power of at least one 
propagator to zero through the use of triangle identities. It is assumed that a tensor reduction 
on the massless 1-loop sub-diagram formed by lines 1 and 3 is performed, such that for A*^^), all 
invariants can be expressed in terms of propagators. Integration-by-parts recurrence relations 
express all integrals of this topology in terms of one primitive integral with non-zero cuts. 
Here and below we do not count integrals that can be calculated by repeated application of 
the closed expressions Eqs.( [A.i| )-( |A.5| ). For the normalization of the primitive integrals we 
define 

^ JlLr'(i-^)r(i + .) ,5,4) 

^ (ttsY r(l-2e) ^ ^ 

where e is related to the space-time dimension as L) = 4 — 2e. 

^It should be mentioned that new approaches to solving integration- by-parts identities are being developed 
as weU M 
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1^(3) (ai,a2, 03,04, 05,06, a?) = i 




d^p d^/fc d^/ 



2]a2 



(5.5) 



Im(I^(3)(l,0,l,l,l,0,l)) =7rs(7V, 



1 / 1 



e \ 2 



+ 



11 

T 



77 15 



+ C2) +e{-— + 
2/ 439 81 69 



C2 + 8C3) 



2 +r:^2 + 50C3 + :^C4)+O(e=^) 



(5.6) 



where for compactness the notation C„ = C(n) has been adopted. 

For 4-loop topology A of Fig.^ one can lower the power of at least one propagator to 
zero through the use of triangle identities. For all 4-loop topologies A-G of Figs.^-^ we 
assume that a tensor reduction on the massless 1-loop subdiagram formed by lines 1 and 3 
is performed, such that one has to deal with only one invariant that cannot be expressed in 
terms of propagators. For topology A this invariant is chosen to he p-Q. Integration-by-parts 
recurrence relations express all integrals of this topology in terms of 9 primitive integrals 
with a non-zero imaginary part. Again we do not count integrals that can be calculated by 
repeated application of the closed expressions Eqs.( |A.i| )-( [A.5 ). 



/a(oi, 02, 03, 04, as, ae, 07, ag, 09, aio, on) 



d^p d^k d^l d^r 




{p-QY 



(r2)"i[(A; + Q)2]"2[(r -k- Q)^]"-^{k^)"i[{k - /)2]"5 



[{k-l- p)2]°6 (/2)a7 (p2)a8 [(^ + p)2 + 2(/ + p) • Q]°9 (P + 21 ■ Q)°10 



(5.7) 
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Im(/A(l,l,l,l,0,l,l,l,l,l,0)) = f(iV,)4 i(-5C4)+O(e0 
Im(/A(l,l,l,0,l,0,l,l,l,0,0)) = Trs{N,) 



+e{-'-W + ¥C2 + ^Cs + f C4) + 0{e' 



Im(7^(l,0,l,l,0,l,0,l,0,l,0)) = ttsHN.Y 



j_^f 3697 119/- 17/- ^ I ^2f 159001 2053 /- 511/- 263/- ^ , /n/'=-3\ 
+^l^~W ~ T2-C2 - -2-Q3 )+£ 36-C2 - -6-C3 - — C4 ) +U{£ ) 



Im(/^(1,1,1,1,0,1,0,1,0,1,0)) = nsiN.y 



+e(^ + fC2-9C3-¥C4)+0(^^ 



Im(/A(l,0,l,0,l,l,l,0,l,0,0)) = TTS^N^Y 



1 1 ^ -L 49 /- \ , 17161 127/- 25/- ^ 
i 1^12 j + 1^18 - C2 j +£(^^32- - T2-C2 - ^C3j 



2/ 260347 2837 



Im(I^(l,l,l,0,l,l,l,0,l,0,0)) 



+^ \^ 648 



36 



A 376 a 319/- \ I 

(,2 3-C,3 r (.4 J + (-^(e 



^ -i]+i(-6) + 



385 I 11 A I ^ 

+ tC2 + C3J 



+,(_651 + 121^^ + 36^3 + ^^^^ + 0(^2) 



Im(/A(l,0,l,l,l,l,0,0,l,l,0)) = 7rs(Ar,/ 



+ ¥C2 + 146C3 + I55C4) + Oie^ 



Im(JA(l,0,l,l,0,l,l,l,l,0,0)) = Trs{N,Y 



1 I 1 

4 



1 / _33 



+ C2) + (-^ + fC2 + 13C3) 



Im(/A(l,l,l,0,0,l,l,l,l,0,0)) = TTsiN.y 



+ f C2 + iFC3 + ^(4) + 

K-i + ^^) + (-^ + TC2 + i3C3 



I ^Z' 4843 I 249/- , 195/- , 125 a A 1 /O/'^-S 
+^(^ 16~ + ~C2 + — C3 + — C4 j + (-'(e 



Topology B has 7 non-trivial primitive integrals with a non-zero imaginary part, not 
counting the integrals that can be considered cases of Topology A. In choosing primitive 
integrals one may prefer in one case to keep a power of (p ■ Q) in the numerator to avoid an 
infrared divergence when the three momenta p, k and I go simultaneously to zero. 
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/^(ai, 02, as, "4, as, uq, aj, as, ag, aio, an) = 

{p ■ Q)"" 



(r2)"i[(A: + Q)2]"2[(r - A: - Q)2]"3(A;2)«4[(A; - p)2]a5 

^ I ] 

[{k-l- p)2]°6 (p2)a7 (/2)a8 [(/ + p)2 + 2(/ + p) • (/2 + 21 • Q)"" i 



(5.8) 
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Im(/B(l,l,l,l,l,l,l,l,l,l,l)) 
Im(/B(l,l,l,l,l,l,0,l,l,l,0)) 

Im(/B(l,0,l,0,l,l,l,0,l,l,0)) 



1 17, 



Trs{N,)^ (-I + C2) +e(-^ + 10C2 + 15C3 



(-^ + 67C2 + 159C3 + IO8C4 )+0{e-' 



Im(/B(l,l,l,l,l,0,0,l,l,l,0)) = 7r{NeY 



F 2 



+ (159-19C2-15C3-^C4) +0(e) 



Im(/B(l,l,l,0,l,0,l,l,l,l,0)) = 7t{N,Y 



1 + (20-C2-4C3) +e(237- 24C: 



Im(/B(l,l,l,l,l,l,0,0,l,2,0)) = f(iVe) 



53C3 - fa) + 0{e^) 

1 / 5 



e \ 4 



C3-3C2ln(2) + 5C3-12C2ln(2) 



-fC4 + 6Li4(l/2) + i 111^(2) ) +0(e) 



Im(/B(l,0,l,l,l,l,0,0,2,2,0)) = ^{N,y 



-3C2 + !C3 + 3C2ln(2) +e -21C2-fC3 



6C2 ln(2) + ^C4 - 6 Li4(l/2) - i ln^(2) ) + 0(e 



Here Li4(x) is the fourth order Polylogarithm [|3T|, Li4(l/2) = 0.517479- • • . 

For topology C one can lower the power of at least one propagator to zero through the 
use of integration-by-parts recurrence relations. For this topology there are three additional 
primitive integrals with a non-zero imaginary part, not counting the integrals that are special 
cases of Topologies A and B. 



Ic{ai,a2, 03, "4, "5, "6, "7, "8, "9, OlO, "ll) 



d^pd^fcd^/d^r 



{kiy 



(r2)"i [{k + Q)2]"2 [{r-k- Q)2]°3 (A;2)a4 [(fe _ p)2] 



"5 



(p2)a6 (I2)a7 (p2 + 2p • Q)"8 [(^ + p)2 + 2{l + p) • Q]"9 (/2 + 2/ • Q) 



"10 



(5.9) 
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Figure 5: Integration topologies D - G. The lines are numbered according to Eqs. (5.10) 
and (Islll) 



- IIC2 + 9C3) + 0{e) 
■ (3C2) + (11 + 26C2 + Cs) + 

fC2-3C3)+6(^-48C2ln(2) 
K3-fC4)+0(e2) 

For topology D one can apply integration-by-parts relations in the form of the massless 
triangle rule to the triangle formed by lines 5, 6 and 7. Further use of integration- by-parts 
expresses integrals of this topology into integrals that can be calculated using the closed 
expressions Eqs. (|A.lD -( |A.5| ) or integrals with no imaginary part, or integrals that are of 
topology A. Therefore no new non-trivial primitives with a non-zero imaginary part appear 
for topology D. 

For topology E one can lower the power of at least one propagator to zero through the use 
of integration-by-parts relations. There is one primitive integral with a non-zero imaginary 
part, not counting the integrals that can be considered as cases of topology A and B. 



Im(lc;(l,0,l,l,l,0,0,l,l,l,0)) = 7:s{N,f 

lm(lc(l,0,l,l,l,l,0,l,l,l,0)) = 7r(iV,)4 ■ 

lm(/c(l,l,l,l,0,0,0,l,l,l,0)) = 7rs(iV,)4 

1/ 3507 _ 

- — (■2 + 
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/^(ai, ^2, "3, "4, "5, ae, ay, as, ag, aio, aii 
d^p dPk dPl d^r 



(p-0 



"11 



(r2)"i [(fc + Q)2]a2 [{r-k- Q)2]"3 (A;2)a4 [(A; - /)2]a5 
1 



[(A: - p)2]"6 (;j2)a7 (p)a8 (p2 + 2;j • (^2 + 2r Q) 



(5.10) 



Im(/£(1,1,1,0,1,1,1,1,1,1,0)) = f(iV,)4 Ji(^_6C4) 

Non-planar topology F has 6 non-trivial primitive integrals with a non-zero imaginary 

part, not counting the integrals that can be considered cases of topologies A-E. In choosing 
primitive integrals one may prefer in one case to keep a power of {p ■ I) in the numerator to 
avoid an infrared divergence when the three momenta p, k and I go simultaneously to zero. 



/p(ai, 02, 03, a4, a5, ae, ay, as, ag, aio, an) 
d^p d^k d^l d^r 



ip-l) 



an 



(r2)"l [{p + k + l + Q)2]"2 [{r-p-k-l- Q)2]«3 (/j2)a4 



[{k + p)2]"5 [(k + ip'^)"^ (^2)"» + 2p • g)"9 (|2 + 2i • Q) 



"10 



(5.11) 



Im(/F(l,l,l,l,l,l,l,l,l,l,l)) = f(iV,)4 

lm(/f(l,l,l,l,l,l,l,0,l,l,0)) = f(Ar,)4 

Im(7^(l,0,l,l,l,l,l,0,l,2,0)) = f(iV,)4 

Im(7F(l,l,l,l,0,l,l,0,l,l,0)) = 7r(Ar,)4 



i(-fC4)+0(e' 
(-fC4)+0(£) 

^(l)+^(l5-C2-C3) 



lm(/f(l,l,l,l,0,0,l,l,l,l,0)) = TTiN.y 



lm(/f(l,l,l,l,l,l,0,0,l,l,0)) = 7r{N,Y 



+ (139 - 19C2 - 17C3 - 9C4) + 0{s) 

i(-l + C2) + (-20 + IIC2 + IIC3) 
+£^-237 + 93C2 + I2IC3 + 37C4) + C(£2) 
i(-2 + C2 + fC3-3C2ln(2))+O(£0) 
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For integrals of topology G one is able to remove enough propagators through the use of 
integration-by-parts relations to make these integrals cases of topologies A-F. Therefore no 
new primitives appear for topology G. 



5.1 The evaluation of primitive integrals 

In this section we will discuss the calculation of primitive integrals on the mass shell using 
a large mass expansion as an intermediate representation. The general theory of Euclidean 
asymptotic expansions that is invoked was developed in Refs. 33 1 and in practice the 
techniques of Ref. were used. In many physical applications the large mass expansion 
is truncated at a finite order, however it has been used in recent works |35,36] also in an 
untruncated form. 

Before treating the more complicated 4-loop integrals we will illustrate some basic features 
of the approach using a simple massive one-loop integral as an example 

I{ai,a2) = j d^k (^2)aip + g)2_^2]a2 = j '^''P [(p + Q)2]«i(p2_M2)-2 ^^'^'^^ 

Again we do not write explicitly +ie in the denominators to keep expressions compact. For a 
large mass, > Qj^g may expand this integral by making an ordinary Taylor expansion 
in the integrand: 

1 1 ,^ /2p-Q + Q2y r(ai +i) 



i—n V 



The remaining integral over p is of the vacuum bubble type and after applying a bubble tensor 
reduction to simplify the powers of (p ■ Q) in the numerator 

rfa2Q2)t n! i;(D/2) 
(p. QY ^ ' 2"(n/2)! r(D/2+n/2) ' ^i^) 

\0 n > odd 

the bubble integral can be performed using Eq. ( [A.l ). In the resulting expression one of the 
two sums can be performed using the summation identity 



^ (-1)'= r(a + A:) _ T{a)T{h-a + n) 



^^k\{n-k)\T{b + k) n\T{n + b)T{h - a) 



(5.15) 
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and obtain the well-known result 

/(ai,a2) = i7r(^/2)(_l)-"i-"2(^2)D/2-ai-a2 



T{D/2-ai) 



r(ai) r(a2) 
^fQ^y T{ai+j)Tia,+a2+j-D/'^l ^ 



j=0 



i,r{^/2)(_l)-"i-"2(M2)^/ 



2— ai — Q' 



T{D/2-ai) r{ai + a2- D/2) 
T{D/2) r(a2) 



X2F1 



ai , ai + a2 — D/2 
D/2 



(5.17) 



where 2F1 is the Gauss hypergeometric function. To obtain the proper value on the mass 
shell Q'^/M'^ = 1 within dimensional regularization, one goes on-shell in an interval of the 
dimension D where the integral is convergent and continues the result in the number of space- 
time dimensions. Indeed, assuming that D is chosen such that the sum is convergent one can 
use the Gauss summation identity 



a , b 
c 



r(c) r(c - g - b) 
r(c-a) r(c-6) 



(C7^ 0,-l,--- , Re(c 



5)>0) 



to obtain the on-shell expression Eq. ( |A.3| ). 

Depending on the values of ai, 02, integral 1(01,02) can have specific IR divergences on 
the mass-shell since [{k + Q)'^ — M^] becomes (A;^ -|- 2Q ■ k). One may use IR power-counting 
in the small k region to determine whether an IR divergence occurs on the mass-shell. Here 
l/{k'^ + 2Q-k) counts for only half the power of an ordinary massless line l/fc^. Since such an 
IR divergence appears only on-shell, it manifests itself as non-convergence of the large mass 
expansion starting exactly on the mass shell. 

Using integration-by-parts recurrence relations it is possible to choose an appropriate 
basis set of primitive integrals for which the large mass expansion is convergent on-shell in 
D = A — 2£. Nevertheless it is instructive to study the behaviour of IR divergent integrals 
near the mass-shell. Let us therefore consider integral 1(1,2). Using the transformation rule 
for the hypergeometric function 



2-t^i 



( a , b 






■) 







r(c) r(c - a - b) 
r(c-a) r(c-6) 



2FI 




. n ^c-a-fe ^(c) T{a + b-c) I 

r(a) r(6) ^^M c-a-b + l 



1 - z (5.18) 



one easily obtains the behaviour near the mass shell 



1(1,2) 



i-K r(l + e) 
2e (M2)i+s 



1 - (1 - x) 



.2, r(i + 2£) r(i-£) 

r(l + e) 



+ 0{l-x) 



(5.19) 
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where x = Q'^ /M"^. On shell the term (1 — x) is nullified within dimensional regularization 
and for x < 1 the threshold term can be expanded in e 



2e 



[1 



il-x) 



-2el 



ln(l -x)+e ln2(l - x) - -e^ ln3(l - x) + O {e^ ln^(l - x)) 

3 



(5.20) 



The logarithms ln*(l — x) in Eq. ( pOD are translated on-shell into a pole. Such logarithmic 
divergences appear for certain individual diagrams of the muon decay calculations but they 
cancel to all orders in e when the external leg corrections are included. One may therefore 
use dimensional regularization and evaluate all integrals at threshold. 

To make the connection with methods where one expands expressions in e at an early 
stage in the calculation, we also expand the T functions of Eq. ( |5.17] ) in e 



r(n + l + e) = n!r(l + e 
where the harmonic sums 5" are defined as 



1 + eSi{n) + ^(S'tin) - S2{n)) + ©(e^] 



n ^ 



(5.21) 



(5.22) 



i=l 



In this way one obtains for integral 1(1,2) 



/(I, 2) 



-ITT r(i + e) 
(M2)i+^ 

-mr(l + e) 

(M2)i+£ 



i=o 



l + 2eSiij) + 2e'Stij) + 



1+j 



ln(l - x) + eln2(l - x) - -e^ ln^(l - x) 



-2C2e^ ln(l -x) + eC2 + 0(£^) + 0{1 - x) 



(5.23) 



in agreement with Eqs.( ^.19 ) and (5.20). However, one needs to consider the divergent loga- 
rithms to all orders in e to obtain the proper on-shell value in dimensional regularization. As 
integration-by-parts recurrence relations are used to reduce the number of different integrals 
that are required, one may conveniently choose a set of primitive integrals that are free from 
threshold singularities. These well-behaved primitive integrals can then be calculated via a 
convergent large mass expansion that can be truncated at a fixed order in e. This is the 
general strategy followed in the more complicated multiloop integrals. The only caveat is 
that, after applying recurrence relations, poles may appear in the coefficients multiplying the 
primitive integrals and these primitive integrals then have to be evaluated to a higher order 
in e than naively expected. 

It follows from Eq. ( 5.18| ) that integral I{ai,a2) has general on-shell divergences of the 
form 



lim 



^(-m-2e) 



♦ 1 L 



(5.24) 
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Figure 6: Integral 1, 1, 0, 1, 0, 1, 1, 1, 1, 0). Thin lines correspond to massless scalar 

propagators and thick lines correspond to massive scalar propagators with mass M. The 
momenta that flow through the propagators are labeled with the line numbers defined for 



main topology B, Eq.(5.8). 



which are nullified^ \n D > 4 + m dimensions. Note that for x > 1 the terms -\-ie in 
the denominators of the propagators that were implicit in Eq. ( ^.12 ) are restored by the 



substitution x ^ x + ie in Eq. ( 5.24 ) to give the contributions from channels that open up 



above threshold, but for the present work this is not needed. 
5.1.1 Example 2 

Let us now consider the calculation of a primitive integral /_b(1, 1, 1, 0, 1, 0, 1, 1, 1, 1, 0) (see 
FigJ^ using the large mass expansion as an intermediate representation. The basic steps 



followed will be similar to the steps in the derivation of Eq. ( 5.17 ). However, the sums 



that the Taylor expansions introduce cannot be reduced at the level of T functions because 



insufficient generalizations of Eq. ( 5.15 ) are known. One is therefore forced to reduce sums 
after expanding the F functions in e. What makes this approach feasible is that the sums over 
harmonic functions, S, encountered in the present calculation can be reduced to expressions 
given in terms of a small number of higher order sums. Eventually one is only interested in the 
value of integrals on-shell. For integrals with more than one massive line one may therefore 
choose to make a large mass expansion in terms of the mass of only one of the propagators, and 
to keep the masses of other propagators on-shell, if this leads to less complicated intermediate 
expressions. 

The problem of reducing harmonic sums also appears in calculations for deep-inelastic 
lepton-nucleon scattering, where one deals with a light cone expansion instead of the large 
mass expansion. Large collections of summation identities of the type that are needed for the 



present work can be found in the literature p7, 38|. Although summation reduction algebra 



■^Studies of the various multiloop integrals that are needed for the muon-decay calculation reveal a slightly 
more general threshold behaviour 

lim^ [(1 



where m is an integer and n is generally the number of independent loop momenta that when they to zero 
give rise to an infrared divergence on-shell. The absence of unregularized singularities (1 — x)~"^ for m > is 
important for the validity of on-shell recurrence relations. 
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is rather involved, relatively compact final expressions for the coefficients of the large mass 
expansion are obtained. On the mass-shell the sums can be reduced to known mathematical 
constants. 

We take the momentum Q off-shell, s = < M^, and perform a large mass expansion 
by expanding the integrand as 



PiPiplplPrPi (pI - M^){pIq - M2) 



2 2 2 2 2 

PiPipiPiPi 



+ 



.i(pi-M2)(pf,-M2) 

1 



pIpIpIpIptpI {pI - M'^)iPw - M"^) 



{P2,P5,P7,Q} 

(5.25) 



{Q} 



where we use the notation that the terms between square brackets are expanded as an infinite 
Taylor series around the point where the momenta between curly brackets are zero. Here 
momentum pi is understood to be the momentum that flows through line i as defined for 
main topology B, Eq. ( |5.8D . After performing the indicated Taylor expansion in Q, the last 
term on the rhs of ( 5.25| ) produces integrals of the 4-loop massive bubble type and therefore 
does not contribute to the imaginary part. The first term on the rhs of ( 5.25| ) produces 
integrals that are a product of a 3-loop massless propagator-type integral and a 1-loop bubble 
integral, all coupled by tensor numerators. After the necessary tensor reductions, one can 
evaluate these integrals in terms of F-functions using Eqs.([01) and ([A1|) . 

A representation for the imaginary part of integral /_b(1, 1, 1, 0, 1, 0, 1, 1, 1, 1, 0) is obtained, 
in this way, as a multiple sum over a product of T functions. After expanding the T functions 



in e as in Eq. (5.21) one can reduce the sums over the functions Si{k) to a basis of independent 
nested sums to obtain the off-shell expression 



s \fc 1 
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(5.26) 



where is defined in Eq. (5.4). On-shell s/M"^ = 1 and the sum over k can be performed 
analytically to give the following result 



Im(lB(l, 1, 1, 0, 1, 0, 1, 1, 1, 1, 0)) = vr {N,f 



1 + (20-C2-4C3) +e(237- 24C: 



53C3-yC4)+0(£') 



(5.27) 
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Figure 7: Integral /^(l, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1). Thin lines correspond to massless scalar 
propagators and thick lines correspond to massive scalar propagators with mass M. The 
momenta that flow through the propagators are labeled with the line numbers defined for 
main topology B, Eq. ( |5.8| ). A factor p-j ■ Q (line 11) in the numerator is not shown explicitly. 



where the various infinite sums that appear are 
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and are known from the literature |3^41|. 
5.2 Example 3 

For integral /_b(1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1) we proceed as in the previous example and make a 
large mass expansion by expanding the integrand 
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(5.33) 

{Q} 



The last term is a 4- loop massive bubble and does not contribute to the imaginary part. After 
performing the indicated Taylor expansions, integrals of the massless propagator and massive 
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bubble types remain that can be evaluated in terms of T-functions using Eqs. ( A.ll )- 
The F- functions are then expanded in e and the resulting sums over the functions Si{k) may 
be reduced to a basis of independent nested sums. Adding the imaginary parts of the various 
contributions in ( 5.33 ) gives the off-shell expression 



Im(lB(l, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1)) = (N,)^ ^ Yl {^) 
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5^« = E(-1)4 



(5.35) 



and for higher order nested sums one defines recursively 

'S'ri rri(j) 



'fc,n,...,m 



r 



k,n,...,m 



r 



(5.36) 



where each of the indices n, ...,m can have a tilde (~) to indicate an alternating sign compo- 
nent. On-shell s/M"^ = 1 and the infinite sum over k can be performed to obtain 



Im(/B(l, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1)) = {N,f f {-^ + O(e0)} 



(5.37) 



To illustrate the delicate cancellation for sjM'^ = 1 between the various divergent sums in 
Eg. ( 5.34 ) let us also quote the result for the sum over k from 1 up to y instead of infinity 
(such that Eq. ( ^.37D corresponds to y ^ oo). 

Im(lB(l, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1)) = (N.f |-C2S2(y) - lc2S^{y) - 52,2(2/) + 3S3,i(y) 
+ ^52,2(y) - 52,i,i(y) + Si{y) [S2,i{y) - 2^3(2/)] + ^5i,i(y) [(2 - S^iy)] } + 0(e°) (5.38) 

Here one encounters in addition to the sums Eqs. (|]2|) - (1^ the known infinite sums 

5^(cx)) = -(l-^)Cfc (A:>1) (5.39) 
S2d^) = TiC4 (5.40) 
^2,1,1(00) = 3C4 (5.41) 

and Si{y) and 6*1,1(1/) diverge logarithmically for large y. 
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5.2.1 Example 4 

As a last example the calculation of integral /^(l, 1,1, 1,1, 1,0, 0,1, 2,0) that produces a 
constant Li4(l/2), where Li^ is the Polylogarithm |^T| is discussed. As before a large mass 
expansion is performed 
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{Q} 



where the last term produces 4-loop massive bubble integrals and does not contribute to the 
imaginary part. After performing the indicated Taylor expansion there remains a product 
of a massless two loop propagator integral and a massive 2-loop bubble integral that can be 
evaluated in terms of F-functions using Eqs.( |A.i[ )-( |A.5| ). After expanding the F-functions in 
£ and reducing the sums over the functions Si{k) to a basis of independent nested sums one 
obtains 
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On-shell the infinite sum over k can be performed to obtain 



Im(lB(l,l,l,l,l,l,0,0,l,2,0)) = f(iVe 



(|C3-3C2ln(2)) + (5C3-12C2ln(2) 



1/5 
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C4 + 6Li4(l/2) + iln4(2) +0(e 



(5.44) 
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for which ah necessary infinite sums are known 
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6 The Theoretical Uncertainty 

With the incorporation of the 2-loop QED corrections in Ag the largest missing theoretical 



contributions come from three possible sources. The first is the hadronic uncertainty |15| 
conservatively estimated to introduce an error of 2 x 10"*^ in Ag which introduces a relative 
error of 10^^ in the extracted value oi Gp- 

By examining the logarithms that appear at tree- and 1-loop levels one would expect the 
leading unknown 2-loop QED corrections to be proportional to 

(a/7r)^(mg/m^) ln^(mg/m^) = 1.4 x 10"*^ which would need a coefficient of roughly 140 to 
introduce a 1 ppm error vciGp- At tree- and 1-loop level the coefficient of the leading logarithm 
is 12. Allowing for a coefficient as large as 24 gives an an estimate of the theoretical error 
from this source of 1.7 X 10"^ in the value oiSGp/GF. 

The 3-loop QED corrections may be estimated in the same way as the 2-loop corrections 
were, before they had been actually calculated, by assuming them to the equal to the known 
leading logarithm of next order. In this case that translates a relative error of 1.4 x 10~^ in 
Gp. 

Overall the theoretical uncertainty arising from missing higher-order corrections should 
not exceed a few parts in 10^. 



7 Experimental Uncertainties 

On the experimental side the accuracy of the known values for the muon lifetime, r^, the muon 
mass and the muon neutrino mass can each exert a significant effect on the extracted value 
of Gp- From Eqs.( |2.5| ) and ( |2.7| ), the change in Gp in response to changes in these quantities 
is given by 

Gf 2t. 2m.^ m?- ^ ' > 
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The measured value of the muon lifetime, = (2.19703 it 0.00004) fis |T[ is currently 
the source of the dominant error. Experiments are planned to reduce this error to ib4ps 
Brookhaven National Laboratory ES] and ib2ps at the Paul Scherrer Institute A new 



measurement is also expected at the Rutherford- Appleton Laboratory |^. It is therefore 
likely that uncertainty on coming from this source will be reduced to somewhere in the 
range 0.5-1 ppm. 

The measured value of expressed in unified atomic mass units is [|[| 

ruf, = (0.113428913 ± 0.000000017) u 

corresponding to an accuracy of 0.15 ppm. A new measurement of Planck's constant, h, 
means that this can now be converted to units of MeV without introducing additional errors. 
The current error on leads to a 0.38 ppm uncertainty on the extracted value ol Gp- Some 
reduction can be expected in the error on the mass whose effect would then be insignificant. 

If m^,^ is assumed to be non-zero then setting it to the current upper bound of m^^ < 
170 keV shifts the extracted value of Gp by 10 ppm. This upper bound is expected to be 
reduced to some where below 30keV by studying the decays of pions in in flight at the 
Brookhaven muon storage ring [47|. This would affect the extracted value of dGp/Gp at the 
level of 0.3 ppm. 



8 Weak Corrections to the Muon lifetime 

The weak corrections to the muon lifetime may be encapsulated in a quantity Ar defined by 
the relation 

^ = -^(l + Ar) (8.1) 

As with Eq.( [4.1| ) the parameters on the left-hand side of Eq.( |8.1| ) represent experimentally- 
measurable quantities and those on the right-hand side are the renormalized parameters in 
whatever renormalization scheme has been chosen. Thus g and are the renormalized 
SU (2) p coupling constant and square of the renormalized W boson mass respectively. 

Ar will be defined through Eq.(^.l|) in such a way that that Eq.(p.4|) is true exactly. 

In an analogous way to Aq, Ar can be expressed as a power series in 



Ar = Ar 



i=0 



where, once again, the index i indicates the power of ar that appears in Ar 



Sirlin |48,|^] has described a strategy that, starting from the full electroweak theory, 
makes the separation of contributions to Aq and Ar automatic at least up to 0{am'j^/M^r). 
In diagrams exhibiting IR divergences, the photon propagator is replaced by 

1 f 1 1 1 1 . , 

+ T^-T^ • (8-3) 
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where it generally convenient to take A = M\y. The term in curly brackets is simply the 
original photon propagator with a Pauli-Villars regulator. It has the same IR behaviour and 
gives contributions that are identical to those of Fermi theory up to 0(am^/A^) and thus 
are contained in Aq. The second term in (|8.3| ) gives contributions that retain the original 
UV behaviour but are free from IR singularities and therefore belong in Ar. It should be 
noted that, in contrast to neutral currents, it is not generally the case that charged current 
processes can be separated into QED and weak parts in a gauge-invariant manner. Sirlin 
has also discussed in general how the strategy may be applied at the 2-loop level. 

Note that the photons in the second term on the right-hand side of ( |8.3D have a mass 
larger than the muon mass such that in contributions to Ar no "photons" of this type can 
appear in the final state. This means that one does not need to perform many body phase 
space integrals for the calculation of Ar and in practice the matrix element can be evaluated 
even at zero external momenta. In contrast the QED corrections to muon decay in the Fermi 
theory do involve many body phase space integrals, but these corrections contain only one 
mass scale, ignoring the electron and neutrino masses, and the relatively simple Fermi contact 
interaction. 

It is known that 

Ar(0) = ^ + of4') (8.4) 

that is due to W propagator effects and shifts the extracted value of Gp by 0.52 ppm. A 
factor (1-|- Ar*^")) is traditionally included in Eq.( |2.4| ) that is used to define Gp and up to now 
this has not been of any significance. It is, however, clearly inappropriate in the definition of 



Gp as it does not arise from the Fermi theory Lagrangian (2.1) and is most properly included 
with the weak corrections. It is striking that at the 1 ppm level experiments performed at low 
energy are sensitive to the finite range effects of the W boson. 

8.1 1-loop Electroweak Corrections to Ar 

The 1-loop corrections to the muon lifetime were first calculated in the Standard Model of 
electroweak interactions by Sirlin |49|. In a general renormalization scheme Ar^^^ can be 



written |51| 



+ -^ o ^— + 2-^ 8.5) 

167r2 24 4 Ml, g ^ > 

^wwi^) 9^ \ J . , 124-7, 1 

+ tI-^ <^ 4 a - In ^ + ^ , In ^ + 6 



^<\2^a^"^ (8 6) 
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where ah parameters are renormaUzed parameters in whatever renormaUzation scheme has 
been chosen. n[y^y(0) is the 1-loop W boson self-energy and sq and cg are, respectively the 



sine and cosine of the weak mixing angle, 6\y, defined as in Ref. |21] so as to diagonalize 
the mass matrix of renormalized Z and photon fields. In Eq.( ^.5p the term in curly brackets 
comes from vertex corrections and the following one from box diagram. Ultraviolet divergences 
cancel in the overall expression. 

8.2 2-loop Electroweak Corrections to Ar 

Certain classes of contributions to Ar are simply related to the corresponding ones in 5p 
because they only enter through the self-energies of the Z and W bosons. In such cases 

Ar = — Up. (8.7) 

Contributions to 6p^'^^ in the limit of large Higgs mass, Mh, were obtained by van der Bij and 
Veltman |52]. Subsequently van der Bij and Hoogeveen |53| calculated corrections to Jp*-^-* 



arising from a heavy fermion doublet. Their results can be used to obtain the well-known 
asymptotic expression 

for large top quark mass, mt, and for which 

xt 



It was this quadratic m^-dependence in Sp that provided some of the strongest constraints 
on the mass of the top quark before it was directly observed. Thus it was the measurement 
of Gf was that allowed the top mass to be successfully predicted from precision electroweak 
data. 

Consoli, Hollik and Jegerlehner showed how to combine ( ^.81) with IPR corrections to 
obtain the asymptotic dependence of M\y on m^. The analogous result for Ar^^^ as defined 
via Eqs.(|3) and is 



At 2-loop order the contributions that are quartic in mt are generated via the Yukawa 
couplings to fermions and therefore come entirely from the scalar sector of the theory. These 



corrections are given, keeping the full Mh dependence, in Refs. |55-59[. 

Since it turns out that rrit is roughly the same order as Mz it is to be expected that the 
0{a'^m^M'^) contributions to Ar^^^ could be similar in size to those that are quartic in mt 
and has been borne out by direct calculation These subleading corrections do not now 
just come from W and Z self-energy diagrams alone and so are not simply related to 6p. The 
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exact sensitivity of Ar^^^ to Mh, without making an expansion in mt has been studied in 



Ref. [|0|. 

None of the corrections mentioned above encounters IR problems. All contribute solely 
to Ar and do not require the invocation of (^]^) in order to separate their QED and weak 
corrections. Recently the 0{Nfa^) corrections to muon decay have been calculated pl[ . 



These are all 2-loop corrections containing a massless fermion loop and (|8.3| ) must be applied. 
In this case, however, the separation is particularly natural. IR divergent contributions come 
from the Feynman diagrams shown in Figs. ^ and ^ 

IR divergent external leg corrections are shown in Fig. ^(a)-(d) with the ' x ' in Figs.^(b) 
and (d) representing the fermionic contribution to the 1-loop photon 2-point counterterm, 
{Q^9iJ.u — QfiQu) 25e^^^^ /e. These can be split into QED and weak contributions using (|8.3|) in 
the manner explicitly described in the appendix of Ref. ||5^ . The QED part is already included 
in Eq.( |2.9D for Aq^'^\ Remarkably for the weak part, the dependence on the separation mass, 
A, cancels in any renormalization scheme between the pairs of diagrams Fig.|Ka) with Fig.|(c) 
and Fig.^(b) with Fig.^d). The combined weak parts of the Feynman diagrams contains a 



simple pole at Z? = 4 that eventually cancels with other divergences from the weak sector [61|. 

The box diagrams, Fig.^, both vanish identically due to a conspiracy in the 7-matrix 
algebra of the fermion currents. The general result is proportional to products of left-handed 
couplings with right-handed couplings; the latter being zero for the W boson. 

In a general renormalization scheme the 0{Nfa'^) box diagrams containing a virtual pho- 
ton and a counterterm insertion on the W propagator produce contributions proportional to 
the IR divergent 1-loop box diagram. 
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^.10) 



and diagrams would need to be treated by applying ( |8.3| ). However in the MS renormalization 
scheme adopted here the sum of the counterterms that appear on the right hand side of 
Eq.(8.10) vanishes. The superscript ^^-^^ denotes the 1-loop contribution to the counterterm 
from a light fermion species [21,61|. 



8.3 Gp in the Analysis of Electroweak Data 

There is an important proviso that must be remembered when using ( ^.3| ) to split up the 
radiative corrections. The resulting sets of QED and weak corrections must be calculated in 
a consistent renormalization scheme. In the present paper, the MS scheme with /i = was 
used for the QED corrections but in other circumstances this need not be the appropriate 
choice. 

In the analysis of electroweak data the numerical values of the renormalized parameters 
of the theory are initially obtained from a set of simultaneous equations for the physical 
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observables as calculated from the full electroweak Lagrangian. To 1-loop |51,^ in the 
Standard Model 



r-' = Yi^|^(l + Ar(°) + Ar«)| (l + AgW+Ag«), (8.12) 
Sp = M| + (5M| + ng^(M|). (8.13) 

where 11^^(0) is the transverse part of the 1-loop Z-7 mixing and Sg'^^^ is the 1-loop U{1) 
coupling constant counterterm. Sp is the position of the pole on the complex plane associated 
with the unstable boson. A form similar to Eg. ( ^.11 ) can be found in Ref. ||5^, Eq.(lla) 



The quantities on the left hand side of Eqs . ( p . 1 1] ) - ( pT3| ) are the physical observables whose 
values are obtained by experiment and those on the right hand side are the renormalized 
parameters and counterterms in the particular renormalization scheme that has been chosen 
which, for the purposes of this discussion, will be assumed to be the MS scheme. The subscript 
has been dropped for the renormalized parameters. For the analysis of electroweak data 
obtained near the resonance it is convenient and natural to take for the 't Hooft mass. 



r 



^ = Mz- Upon solving the equations Eqs.(8.11)-(8.13) one finds Or ~ 1/128 and thus the 



running of a is automatically incorporated in a self-consistent manner at tree-level which 
avoids the need to resum large logarithms as is required when the on-shell renormalization 
scheme is used. It also has the consequence that the quantity in curly brackets in Eq.(8.12) 



must also be evaluated at /x = Mz- This can be implemented by substituting = ae(Mj) 
in the formulas ( |2.8| ) and ( p.9| ) for Aq^^^ and Ag(^). The net result is that the quantity 
in curly brackets effectively defines a running Fermi coupling constant, GF{^^) for which 
Gf{Mz) = 1.16639 X 10~^ GeV~^ very close to quoted value obtained for jj, = but 
including only 1-loop QED corrections. 
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A Elementary Multiloop Integrals 

In this appendix expressions are given for a number of integrals in dimensional regularization 
for which exact analytic results are known. 

For a 1-loop massive bubble integral one has 



/ 



ijp' — w? + ie)"i {p^ + ie] 



02 



^ ' ^ ' T{ai)T{D/2) ^ ' 

A 1-loop massless propagator-type integral has the simple form 



•l-D (g2-)g/2-ai-a2 ^"(-P/^ ~ «l)r(I^/2 - Q2)r(Ql + 0^2 " -D/2) 2) 

r(ai)r(a2)r(Z) - ai - 02) 
and a compact expression for this integral with a general tensor numerator can be found in 



Ref. [63|. For a 1-loop on-shell propagator- type integral one can obtain 
dPp 
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and for a 2-loop bubble integral with one massless and two massive lines one finds [p2, 54 

d^p d^k 



(p2 - m2 + ie)"i [{p + kY + ie]°2 [k"^ - + ze)"3 

r(ai)r(a3) 
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also for this integral with a general tensor numerator a compact expression is known |65|. 
Several more simple cases follow by using these expressions recursively as the powers 01,02 
and 03, are allowed to be non-integer, possibly containing D. In this way one can obtain for 
instance 

d^p d^A; 



(p2 + ie)ai (A;2 - m? + ze)°2 (A;2 + ie)"* [{p + k^ + ie]"3 
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(A.5) 



37 



Ai(2) Al(3) A2(3) A3(3) 



Figure 10: Two and 3-loop diagrams whose cuts give contributions to the muon decay rate. 



B Results for Individual Diagrams 

The results for the 4-loop diagrams Al - G2 (see Fig.|), the 2- and 3-loop dia grams 

Al(2) 

- A3(3) (see Fig. 0) and the mass counterterm diagrams XI, X2 (see Fig. |l^) are listed 
below. The diagrams are calculated in a general covariant gauge for the photon field, i.e. 
the photon propag ator reads -i [g^"" - Cq^q^/iq^ + ie)] /{q^ + ie). The average over muon 
spins is taken in the standard way by closing the 7 string of the external muon line with 
(^ + m^)/2. Diagrams are calculated in the Fierzed form (charge retention order) with 
only the vector part of the V— A interaction; restoring the axial vector contributions yields 
an overall factor 4 that is not included in the diagrams below. Furthermore, diagrams 
A1(3),A1,A2,A3,A6,A7,B1,C1,C2,D1,D2,G1,X1 and X2 still have to be multiplied by a sym- 
metry factor 2. The regularization parameter e is defined through = 4 — 2e with D the 
space-time dimension. 

For the imaginary part of the 2 and 3-loop diagrams in Fig. |l^ one has 



Im(Al(2))/N(2) = l + e(^f^+e^(^^-fC2^+Oie^ 
Im(Al(='))/N(3) = i(|-|e) + f^-8C2-^e) 



+, I 10855 _ 218^^ _ 24^3 - 2122^ ^ ) + 0{S' 



Im(A2(3))/N(3) = 1(^1 + + (275 ^112^ 



+e(^-24C2 + ^e-12C2e)+0(£^ 



Im(A3(3))/N(3) = + + + 

+ J _ 2122 ^ I2C2 + 12C2e ) + 0{e' 
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here and below the normahzation factor for A;-loop diagrams is 

^(k) ^ Gls^ /a,xfc-2 / (47rr /^^^ T\l - e)T{l + e) \'' 
10247r3 Utt/ V r(l - 2e) ) 

with s = m^. For the 4- loop diagrams in Fig || one obtains up to 0{£) 

im(Ai)/NW = ^(|-|e + r) + i(!f-8C2-iFe + 8C2C + ^e^) + (^ 

-iff Cs + 32C4 - if^e + '-fC2^ + f CsC + ifi^e^ - f C2e2) 
im(A2)/NW = ^(-iF-fe-ie^) + ^(-^-88C2-if^e-8C2e-^e^) 

,/ 382951 1124 ^ 728 a 98471 t , 256 a t 88 a t 15745 e2 , 
+ I 2T6 3-'^2 - — Cs - Tor 5 + — C2§ - T'^S? - ^16-? + 

im(A3)/NW = + + + + + 

I / 340409 1136 /- , 290629 1 608 /- t , 15745 12 20 a t2\ 
+ (^08 3-'^2 + ^T^? - — C2e + TTT? -T^24 I 

Im(A4)/NW = + + + + 

+(-^ + ¥C2 + - f C2C + - f C2e^) 
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Im(A5)/N(^) = ^(i + i^ + ie^)+i(f^ + ^e + m^) 

I / 59759 68/- , 40295 1 , 64 a t , 20689 ^2 32 /■ ^2] 
Im(A6)/N(^) = ^(|-|^-|^2^+l(379 _i6^^_197^_8^^^_l|5^2^ 

+ - 226C2 - ^C3 - i^re - 68C24 - if CaC - ^ST^^ + f (2^^) 

im(A7)/NW = ^(-| + |e-|e2) + i(-m+8C2 + ^e-8C2e-We^) 

,/ 25177 I nnA j_ 88 /- , 20479 1 302 /- c 88 a t 15781 c2 , 32 a ^2 \ 
Trv.f'RI V1\t(4) _ 1/ 8 ,4/:i lt2\ I / 73841 4808 a i 1024 ^ }„ir,\ , 1700 a 56a 

Im(Bl)/N^ - 7(-3 + 3^+3?J + l- T86 2rC2 + — C2 ln(2) + -^(3 - -^(4 

+W^+fc2^+we) 

Im(Cl)/N(^) = ^(i-|C + r)+i(^ + 8C2C-^^ + We^-8C2) + (-™-fC2 

-f C2 in(2) + 6OC3 ^C2? + f C3e + we - fc2e) 
im(C2)/N(^) = ^{f + ic-ie) + i[^ + f^-we) 

I I 159511 112 a loa/- 8507a 16a a 15745 a2 , 20 a a2 \ 
+ 1 ~324 rC2 - 1^8C3 - 308"^ - -3-C2^ - ^le-? + tC2? I 

Im(C3)/NW = i(-| + fC2) + (™ + ^C2-if^C3) 
Im(C4)/NW = ^(-f)+i(-iF) + (-^ + ifC2) 
Im(C5)/N(^) = i(|) + (ff) 

im(Di)/NW = ^(^i-ic+ie) + i[^-8C2-'-^^+8C2c+we) 

+{w- ^C2 - f cs + f^c4 - + TC2C + ifcae + iiTc^ - fC2e) 

Im(D2)/N(^) = ^(^-| + |e-|e2)+i(-lM + 8C2 + ^C-8C2e-ire^) 

+ ( - ^ + ¥C2 + ^C3 + - i22C2e - i§^c3C - ifir + fC2e) 



40 



im(D3)/N(^) = + + 



Im(D4)/N(4) 

Im(D5)/N(^) 

Im(D6)/N(^) 
Im(D7)/N(4) 
Im(D8)/N(4) 
Im(El)/N(^) 

Im(Fl)/N(^) 
Im(Gl)/N(^) 
Im(G2)/N(^) 



+ 



1 / 2 



14887 I 64 
108 



I D4/- I 8078 1 

+ — C2 + —4 



128 



C2i 



4039 12 I 64 



1/161 m 

9 



s \ 18 



+ 



/ 16819 _ 32 



C2 



16819 , 



i + f C2e + 



18 



16819 c2 32 



\^ 216 3 "'^ 108 ^ 3 ^ 216 3 



1 / 4 



* C _i_ 4 t2 \ I 1 / 155 310 t I 155 e2 
i^l3-3? + 3? e 



+ 



/ 15637 _ 64 A _ 15637 t , 128 

I 108 3 ^2 54 ? 3 '^^'^ ' 108 



8 _L 32 a \ I / 1585 _i_ 3088 a , 544/- 
3 + + - — + — C2 + — C3 



3/+t 



410 



+ 



22937 I 16 



54 



, IDA 
+ -C2 



^(|-|^ + r)+i( + ^-16C2-^? + 16C2C + W) 



+ 



^,2 9 ^3 3 1,4 57 s ' 3 s2<; 3 ',3<; ins s ^ 152^ 



448, 
3 



9389 , 
27 



604 
3 



208, 
3 



15637 c2 



108 



( 68945 _ 16700, 
I 324 

/ 38200 16142/- 1024 a 680 a -inA/- , 56 a , 232 a a oqa a , 

I ^43 8]~'^2 + — C2 in(2j — g-(;3 - 1U4(;4 + -9-? + —<~,2t, - O^C3? + 

- i - iff C2 + ^C3 - - f C2e + 32C3e - le + 
i(-f-fe-|^^) + ("^ + fC2"64C3 + fe-64C2e-!fe^) 



4rf 
3? 



For the mass counterterm diagrams XI, X2 of Fig. 11 we obtain up to 0{e) 



Im(Xl)/NW = ^(^40 + 20e^ +i(^ + 96C2 + ^e^ 



^ - 196C2C + ™e + 416C2 + 288C3 



Im(X2)/NW = ^(^- 52 - 264)+i(-if^-if^e) + 



i^ + 416C2-^e + 208C2? 
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^ ^ (5m^ ^ (5m^ 

XI X2^'^'^'^'^^ 



Figure 11: Diagrams XI, X2 that have a muon mass counterterm insertion bm^ given in Eq. 



Finally, for the external muon-leg renormalization we need to add IPR diagrams that are 

(B.2) 



generated by multiplying the IPI diagrams of Fig. 10 by the gauge independent factor 

1 



1 + Z(a,e)+Z^(a,e) + 



1 - Z(a,e, 

Here Z{ol^ e) follows from considering the proper muon self-energy 

^{m^, q) = {(/l- m^jSi(m^, q^) + m^S2(m^, q^) 

and calculating 



d 

Z = Si(m^,m^) + 2m^^S2(m^,g^) 



(B.3) 



(B.4) 



After evaluating the necessary diagrams contributing to the self-energy and the derivative 
one obtains 



Z(a,e) 



- + 4 + e(8 + 3C2) + 0(e^) 



where 



1 /45 19 \ /79 

-- f — - 9n/ - —uh] - ( Y + 87C2 - 96C2 ln(2) + 24(3 
-yn/ - 12C2n/ - ^^riH + 24C2n/^^ + 0{e) + 0{a^) 

^ (ar\ (47r)" r^(l-£)r(l + £) 



v47r 



\ra. 



r(l - 2e) 



(B.5) 



(B.6) 



and nj labels the number of light fermions (in the present case only the electron, n/ = 1) 
and uh labels the number of fermions with a mass equal to the muon-mass (in the present 
case only the muon itself, nn = 1). 
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C The 0{a) Electron Spectrum 



The complete 0{a) corrections to the electron energy spectrum in muon decay were first 
calculated by Behrends et al. Q who gave the results in terms of hyperbolic functions. As 
was pointed out by Berman |^6| this calculation did not take into account longitudinal degrees 
of freedom for the photon that must be included when infrared regularization is done using a 
photon mass. The necessary correction term was provided in Appendix C of Ref. ||Tl|. The 
calculation was repeated and confirmed by Grotch but dropping terms of 0{am'l/rin?^) 
in well-defined and reconstructible places. A related calculation, in the context of QCD 
corrections to heavy quark decays, has been performed in Ref. The expressions given 

were used by Nir ||l^ to obtain the mass-dependent corrections given in Eq. (|2.8| ). 

We have taken the results of Ref. Q and rewritten them as function oi z = y/xe^ where 
X = m1/m?^ and 6 is the variable, introduced in Ref. |^, defined by cosh.0 = pe • p^/(^^^'e"^/i)• 
In the rest frame of the muon, cosh 9 = Ee/me and hence = (Se +l'e)/"ie where and pe 
are the energy and 3-momentum of the electron respectively. The electron energy spectrum 
is then 



P{z)dz = ^^[Po{z) + ^P^{z)) dz (C.l) 



in which 



2(2^ — x)^ 

Po{z) = ^ ' (2g;^ - 3x(l + x)z + %xz^ - 3(1 + x)z^ + 2z^) , (C.2) 

p^i^z) = _^i£__^ (5^2 _ I0x(l + + (11 + 8x + llx^)^^ - 10(1 + -F5z^) 



i^-xf 

2 



2z5 



z 



(4x^ - 3x(l + 3x)z + IQxz^ - 3(1 + 3x)z^ + Az^) In: 



+ - — r^(4x^ - 3x^(5 - 3x)z + 12x(l + 4x - 2x^^-2 



z- 



Qz^ 

+{5 - 87x + 9x^ + 5x^)z^ 

+12(1 + 4x - 2x2)z4 - 3(5 - Zx)z^ + 4^^) In — (C.3) 

X 



z5 



(2x'' - 3x(l + x)z + 8xz^ - 3(1 + x)z-^ + 2z^ 

+ In (l-^A- + ln(l 
z^+x V zJ z^+x 

, 1 - z , 2; , Z - X 

+ m In — h m z In 



Z X z^ 



-2U2 (1 - J) + 2Li2 (1 - ^) - 2Li2(l - 

When integrated in the interval z G (y^, 1), Po{z) yields the expression for Aq^^\ with 
y = ml /m? = 0, as given in Eq.(p.7[) and (a/vr) Pi{z) gives Ag'^^) of Eq.(^^). 
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The 1-loop QED corrections to the differential decay rate of radiative muon decay, fi 
i^fi'^el, are given in Ref. |69|. 



D The Branching Ratio for fi~ 



(2) 

The contribution Af/^jj^ to the muon inverse Ufetime, given in Eq. (|2.13D , includes pieces that 



come from the processes ^~ — > e~v^i'e, e~f^Pe7 and e~f^Pee e~. In order to obtain the 
branching ratio for the last of these processes alone, the contributions from the other two 
need to be subtracted. Contributions from the process ^~ — > e~f^^'e7 can be eliminated 
by adopting the on-shell renormalization scheme since fermion loops on real photon lines, 
with = 0, are forced to vanish. By the same argument used in section the conversion 
from the MS renormalization scheme (/i = m^) to on-shell scheme can be performed for the 



contribution of electron loops ( 2.13 ) by adding a term 



=(1) 
'OS 



5e.. 



ill 

'MS 



a\ 2 
vr. 



--C(2)l In^ 
24 ' m 



(D.l) 



The contribution in Eq.(2.13) coming from the process e^v^Ue, in which there is 

virtual photon line containing an electron loop, may be calculated using the dispersion relation 
methods of Ref. [|l5| . Note that the use of dispersion relations naturally invokes a subtraction 
at = and is therefore consistent with the on-shell renormalization scheme. ^From Eq.(9) 



of Ref. [15 1 the contribution from Feynman diagrams with virtual photons containing electron 
loops is given by 



Ar 



a 



had 



dz 



3vr ./4p z 



R{mlz)^T{z) 



(D.2) 



where Ar(2;) is given in Eq.(7) of that work and in the present case p = ml /mi and 



R{mlz 



1 + 



2p 



1 



4p 

z 



(D.3) 



The integration may be conveniently divided into two regions, z £ [4/3,4] and z G [4, oo]. In 
the latter interval the integral contains no singularities in the limit me and p can be 
safely set to zero. Making the substitution u = \J\ — A/ z, the integral in this interval is 



AT 



[4,oo] 



K{u)du 



(D,4) 



Vvr/ \5184 36^^ V216 -^^"^ ' ' ^ ^ ' 



-To ( ^) 0.0421308 



C(3) 
(D.5) 
(D.6) 
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where K{u) is given in Eq.(ll) of Ref. |^]. On the interval z £ [4/?, 4], the integral does 
exhibit singularities as — > 0. Setting u = iv\J — 1 in the integral over the first interval 
gives 

= -Fo (^)%.47056 (D.8) 
where the integration has been performed numerically. As mentioned above, the integral 



( D.7 ) contains terms that behave like In"/? for non-negative integer n < 3 that are singular 
in the limit nie — > 0. It will also contain terms that vanish in that limit but such terms were 
discarded in the original calculation of the kernel K{u). Thus the above result is correct only 
in the terms that do not vanish as nie — > 0. 

Adding Eg. (p. 51) to Eg. ( 2. 13 ) to convert to the on-shell scheme and subtracting ( p.8| ) to 
leave only the contribution of p^ e~v^U(,e^e~ gives the branching ratio 



..x n 25361 _ 25 ^ ^ ,^ _ 2 A 2 

Vtt; \ 5184 9 ^^^ V2I6 3^^ 7 9 ^ ^ 



_ /a\2 f 25361 25 



->) + (|-C(2))ln|-9.47056 ^ (D. 



a 



„ 2 

-) 6.30028 (D.IO) 

TT/ 

= 3.40 X 10~^ (D.ll) 
in good agreement with the experimentally measured value of (3.4 it 0.4) x 10~^ Q. 
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